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1 Introduction 

Recently, hydrogen is paid much attention as new primary energy source because of the 
depletion of fossil fuels and environmental issues such as global warming. There are three 
basic research challenges - production, storage, application - for a "hydrogen economy" . 
Our study in this paper is related to hydrogen storage. As is emphasized in the recent 
reports on basic research challenges for hydrogen storage, its high efficiency is a key factor 
in enabling the success of the hydrogen economy (l||2J. Hydrogen storage systems must 
exhibit following properties: appropriate thermodynamics, fast kinetics, high storage ca- 
pacity (more than 10 wt %), effective heat transfers, high volumetric densities, long cycle 
lifetime, safety under normal use |3|. To implement these properties, solid-state storage 
is useful. Metal hydrides, chemical storage materials, nanostructured materials are well 
known for effective solid-state storage systems. 

Among them, we investigate aluminum hydrides in the form of aluminum clusters. 
More specifically, we deal with an aluminum tetrahedral cage (AI4) and its hydrides which 
were recently observed in experiment and confirmed to have enhanced stability plIHl. Al- 
though the clusters found in Refs. ^[5] are produced by vaporizing aluminum metal in 
hydrogen gas and do not lead to immediate application for hydrogen storage material at 
this stage, it would be useful to study their properties theoretically to further explore pos- 
sible connection with more realistic materials for hydrogen storage. It is also interesting 
to study the tetrahedral cage in the aspect that it is the fundamental structure of an alu- 
minum icosahedral cage (AI12) [6]. In this paper, we report on the structures of aluminum 
hydrides which can be constructed by adding H2 molecules to an AI4 tetrahedral cluster. 

We also evaluate and discuss the chemical properties of aluminum clusters and hydro- 
gen adsorption by using a novel method of the electronic stress tensor based on the Regional 
Density Functional Theory (RDFT) and Rigged Quantum Electrodynamics (RQED) [?]- 

[19). 

This paper is organized as follows. In the next section, we briefiy explain our quan- 
tum chemical computation method. We also describe our analysis method based on the 
RDFT and the RQED, and in particular we define our bond orders and regional chemical 
potential. In Sec. [3| we discuss our results. Sec. |3.1| shows structures of hydrogenated AI4 



tetrahedral clusters and discuss their binding energies. In Sec. |3.2| we analyze the struc- 
tures using the electronic stress tensor and our bond orders. We summarize our paper in 
Seed 

2 Theory and calculation methods 

2.1 Ab initio electronic structure calculation 

We perform ab initio quantum chemical calculation for several clusters of aluminum hy- 
drides derived from an AI4 tetrahedral cage. In this work, calculations are performed by 
Gaussian03 program package [20] using density functional theory (DFT) with Perdew- 
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Wang 1991 exchange and correlation function (PW91) f2Tp2]. The spht-valence triple-zeta 



6-31 IG** basis set 23-251 with polarization functions has been used. Optimization was 



performed without imposing symmetry. The visualization is done using Visual Molecular 
Dynamics (VMD) |26| and PyMOL Molecular Viewer programs |27|. 



2.2 RDFT analysis 

In the following section, we use quantities derived from the electronic stress tensor to ana- 
lyze electronic structures of hydrogenated AI4 clusters. This method based on RDFT and 
RQED [7-13,19 include useful quantities to investigate chemical bonding and reactivity 
such as new definition of bond order and regional chemical potential [l4- 16 . We briefly 
describe them below. (For other studies of quantum systems with the stress tensor in a 
slightly different context, see Refs. [28] - [39] .) 

The basic quantity in this analysis is the electronic stress tensor density r ^{f) whose 
components are given by 
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where {fc,/} = {1,2, 3}, m is the electron mass, '0i(r) is the ith natural orbital and Vi is 
its occupation number. 

Taking a trace of r ^{f) can deflne energy density of the quantum system at each 
point in space. The energy density e^{r) is given by 



4(r-) = ^E-''' 



(r). 



(2) 
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We note that, by using the virial theorem, integration of e^{r) over whole space gives usual 
total energy E of the system: J £^{f)df= E. 

Regional chemical potential fiR \7\ is calculated approximately using sf^f) [14]. 
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where n(r) is the ordinary electron density at r. Since electrons tend to move from high 
fiR region to low jur region, the distribution of fiR maps the chemical reactivity. 

Now, we deflne bond orders as £^{f) or /llr at "Lagrange point" 14 . The Lagrange 



point Tl is the point where the tension density r (r) given by the divergence of the stress 



tensor 
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vanishes. Namely, r^^{rL) = 0. r^{f) is the expectation value of the tension density 

operator T^{f)^ which cancels the Lorentz force density operator L{f) in the equation of 
motion for stationary state [ll]. Therefore, we see that T^{f) expresses purely quantum 
mechanical effect and it has been proposed that this stationary point might characterize 
chemical bonding |14|. Then, our newly defined bond orders are 



^rwi^L)' 



and 



One should note normalization by the respective values of a H2 molecule calculated at the 
same level of theory (including method and basis set). 

We use Molecular Regional DFT (MRDFT) package |40j to compute these quantities 
introduced in this section. 



3 Results and discussion 

3.1 Structures and stability 

The bare tetrahedral optimized structure of AI4 is shown in Fig. [l] (a). We found that all 
the six Al — Al bonds have an equal length of 2.74 A to a great accuracy. AI4 is considered to 
have a structure very close to a regular tetrahedron. We note that this regular tetrahedral 
structure is stable only for high spin state (multiplicity = 5). Also it should be noted that 



the global minimum of AI4 cluster is planar rhombus with multiplicity 3 41,42 , which 
has lower energy by 0.5 eV. Below, we investigate the structures when hydrogens adsorbs 
to this tetrahedral structure. 

We first considered the adsorption of two hydrogen atoms. We examined many com- 
binations of adsorption sites and multiplicities 1, 3 and 5. We found that AI4H2 with two 



hydrogen atoms at terminal sites as shown in Fig. [T] (b) and with multiphcity 3, has the 
lowest energy. 

We further added hydrogen atoms to this structure. Fig. [l] (c)-(g) show structures of 
most stable isomer of Al4H^ (n = 4, 6, 8, 10 and 12) we have obtained. We have tried to 
adsorb more hydrogen, without success. We could not find a stable structure for AI4H14. 
Thus, we believe AI4H12 is the saturated structure of the AI4 cluster. 

Here, we comment on the comparison with the structures which were reported in 
literatures. The structure of AI4H4, Fig. [l] (c), is consistent with Ref. (5,43 and A^Hg, 
Fig. [1] (d), is consistent with Ref. [5|. As for AI4H12, there are several literatures which 



investigated stable structures for this as a tetramer of alane AIH3 143-46 . This structure 
also attracts interest because of its high hydrogen storage capacity of 10.0 wt%, which 
exceeds the target value of a hydrogen storage system specified by some technical report [3] . 
Ref. [43] has pointed out that structures in Refs. |44j and j45j are not stable in light of 
recent quantum chemical computation and reported a structure with S4 symmetry as the 
global minimum. Our results. Fig. [l] (g), agrees with the structure found in Ref. (43) . 
Ref. (43] also reported stable structures of AI4H10, including the structure we have found 
as Fig. [ij (f ) . They have shown that there is a structure with lower energy but since that 
structure is chain-like, we adopt the structure of Fig. [l] (f ) as the one which is derived by 
adding hydrogen atoms to Fig. [l] (e). 

In order to investigate the stability of these structures, we here define two types of 
binding energies (B.E.). The total B.E. of Al4H^ is 

AEtotai = E{Amj - [4E{Al) + mE{li)], (7) 

where E(X.) is the energy of X. AE'totai represents the sum of the strength of all the bonds 
existing in the molecule. The average B.E. of H atoms is defined as 

AEn = - {^(AUH^) - [^(AU) + mE{E)]} , (8) 

m 

which represents the strength of Al — H bond per one hydrogen atom. We use the structure 
of Fig. [l] (a) to calculate E^AU). 

The results for these two types of binding energy are summarized in Table [T] and 
Fig. [2J As is shown in Fig. [2] (a), AE'totai decreases as more hydrogens are adsorbed. This 
indicates that adsorption of hydrogen stabilizes the cluster. On the other hand, as we plot 
in Fig. [2] (b), AEa is almost same for AI4H2 and AI4H4 but increases as more hydrogens 
are adsorbed. This is considered to be due to the increased hydrogens at bridge sites for 
larger clusters than AI4H6. It should be stressed that lower total B.E. and higher average 
B.E. of H atoms are favorable for hydrogen storage systems. This means that hydrogen 
atoms turn into state such that they are easily desorbed, while clusters become more stable 
as hydrogen are adsorbed. 



3.2 Stress tensor analysis of chemical bond 

In the previous subsection, we have shown structures of aluminum hydrides derived from 
the AI4 tetrahedral cage and discussed the stabihty of these clusters by the usual binding 
energy. In this section, we discuss their chemical bonds using RDFT analysis introduced 
in Sec. 12. 2[ 

Actually, we have already used the RDFT concept to draw Fig. [l} There, we draw 
bond lines when the Lagrange point is found between two atoms. The Lagrange point 



is the point at which tension density vanishes (Sec. 2.2) and considered to be suitable to 



define chemical bond [T4]-[T^. Two types of bond orders, b^ (Eq. ([5)) and 6^ (Eq. ([6 
are computed and summarized as functions of bond distance in Fig. 3) The b^ describes 
bond strength in relation to bond in H2 molecule and how particular bond contributes to 
lowering of total energy of a structure, whereas 6^ shows bond electrophilicity in relation 
to H2 molecule as reference bond. 

In the figure, we distinguish diflFerent types of bonds. Al-Al bonds with and without 
bridging hydrogen are denoted by "Al-Al(b)" and "Al-Al". Al-H bonds at terminal sites 
are denoted as "Al-H(t)" and those at bridge sites are denoted as "Am(b)". We find 
good correlation between bond distance and our bond orders as has been found in other 



molecules |[14-16 . b^ and 6^ basically exhibit similar features regarding the correlation 
between bond length in the sense that the slope for Al-Al bonds are larger than Al-H 
bonds. There is a subtle diflFerence in AHI bonds between b^ and 6^. For 6^, Al-H(t) and 
Al-H(b) together make a single slope but there seems to appear three families of slopes 
for 6^. Since 6^ concerns the regional chemical potential, the difference may refiect the 
chemical reactivity of these bonds. 

We can further examine features in chemical bonds by analyzing the electronic stress 
tensor (Eq. ([l])). For example, hydrogen bridged Al-Al bonds are investigated for AI4H8, 
AI4H10 and AI4H12. Figs. [4], [5] and |6] plot the largest eigenvalue of the stress tensor and 
corresponding eigenvector on a plane which includes three atoms which participate in the 
bridging bonds, respectively for AI4H8, AI4H10 and AI4H12. Since the eigenvectors have 
three spatial components, namely they are 3D objects, we express them by projecting on 
the plane. In all these clusters, there are Lagrange points between Al and H. There is a 
region with positive eigenvalue of the stress tensor (tensile stress) between them, which is 



typical for covalent bond involving H atom 12, 18 . Also, there is a fiow of corresponding 



eigenvectors connecting Al and H which indicates a formation of strong bonding. 

For AI4H8, which has a Lagrange point between Al atoms, also has a fiow of eigenvectors 
between them. In contrast to the case of Al-H bond, the eigenvalue of the stress tensor 
turns out to be negative (compressive stress) indicating a diflFerent nature of bonding. In 
the case of AI4H10, we did not find a Lagrange point between Al atoms. In Fig. [5} there is 
a similar fiow of eigenvectors as in Fig. [4] between Al atoms, but this structure is shifted 
ofi[ the region between Al atoms likely to be because the H atom came close to Al atoms. 
Therefore we may conclude there is not so much direct interaction between Al atoms as to 
call bonding. We can reach similar conclusion for the case of AI4H12. Fig. [6] shows similar 



structure but it is further away from the region between Al atoms and in fact (when we 
see the flow in 3D) the flow is connected to H(ll) atom which locates at shghtly off this 
plane. This indicates that Al-Al bonding is completely disrupted by the existence of H 
atoms. 

Another example of the electronic stress tensor analysis is provided for AI4H10 in Fig.[7| 
This concerns the somewhat radical structural change between n = 0,2,4,6,8 and n = 
10, 12 of A^H^. As shown in Fig. [l} in terms of Lagrange points, four Al atoms form a 
tetrahedral cage for n < 8 and but the cage seems to be broken for n > 10. This is also 
seen in a jump in the average Al-Al distance (Table [l]). We can conflrm this by analyzing 
the stress tensor density. As is shown in Fig. W\ there is a region between Al(2) and Al(3) 
in which eigenvectors go perpendicular to the plane (so that expressed by dots), showing 
total disconnection of these atoms. 



4 Summary 

In this paper, we investigated the structures of aluminum hydrides derived from a tetra- 
hedral aluminum (AI4) cluster using ab initio quantum chemical calculation. We reported 
stable structures of A^H^ (n = 0,2,4,6,8, 10 and 12), which include structures already 
found in the literature who had investigated the hydrogenated aluminum clusters from 
other aspects. We calculated binding energies of the aluminum hydrides and found inter- 
esting properties as hydrogen storage material: stability of the clusters increases as more 
hydrogen atoms are adsorbed, while stability of Al-H bonds decreases. 

We also analyzed and discussed the chemical bonds of those clusters by using the 
electronic stress tensor. The bond orders deflned from energy density and regional chemical 
potential (which are in turn calculated from the stress tensor) are shown to have a good 
correlation with respect to the bond distance and to be able to distinguish types of bonding 
to some extent. As far as metallic elements are concerned, our bond order analysis had 
been only applied to Pt clusters |16| so the present analysis can be useful basis for further 
research using our stress tensor based analysis. This is also true for the eigenvalue and 
eigenvector analysis of the stress tensor. We have found that Al-H bonds have a positive 
eigenvalue (tensile stress) at the region between the atoms where as Al-Al bonds have 
a negative value. This indicates that the stress tensor can be a powerful tool to classify 
chemical bonding and may provide a deeper insight into the nature of chemical bonds. 
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Table 1: Total B.E. AE'totab average B.E. of H atoms A£'h and mean nearest-neighbor 
bond lengths rfx-x of Al4H^ (n = 2, 4, 6, 8, 10 and 12). H(t) and H(b) denote a hydrogen at 
the terminal site and the bridge site respectively. 





AEtot.i (eV) 


A^H (eV) 


dAl-Al (A) 


C?A1-H(t) (A) 


C?A1-H(b) (A) 


A14 


-5.59 


— 


2.74 








A14H2 


-11.5 


-2.94 


2.66 


1.61 





A14H4 


-17.3 


-2.94 


2.60 


1.61 





A14H6 


-22.7 


-2.84 


2.63 


1.60 


1.74 


A14H8 


-27.6 


-2.75 


2.69 


1.59 


1.77 


A14H10 


-32.2 


-2.66 


3.11 


1.59 


1.74 


A14H12 


-36.6 


-2.59 


3.47 


1.59 


1.72 



(a) 




Figure 1: Optimized structures for AI4 tetrahedral cage and its hydrides A^H^ (n = 
2, 4, 6, 8, 10, 12). The bonds are drawn at which Lagrange points are found and our energy 
density based bond order b^ is shown. 
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Figure 2: Total B.E. AE'totai (panel(a)) and average B.E. of H atoms AE'h (panel (b)) of 
AUH^ (n = 0, 2, 4, 6, 8, 10 and 12). 
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Figure 3: Energy density bond order h^ (panel (a)) and chemical potential bond order 
hy, (panel (b)) as functions of bond length. Data are taken from the structures of A^H^ 
[n — 0, 2, 4, 6, 8, 10 and 12) as shown in Fig. [!} Al-Al bonds with bridging hydrogen are 
plotted with red crosses and those without are plotted with magenta plus marks. Al-H 
bonds at terminal sites are plotted with blue asterisks and those at bridge sites are plotted 
with green squares. 
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AI4H8 





Figure 4: The largest eigenvalue of the stress tensor and corresponding eigenvector of 
AI4H8 on the right panel. They are shown on a plane which includes three labeled atoms. 
As for the eigenvectors, the projection on this plane is plotted. The positions of these 
atoms are shown by the circle dots. Parenthesized numbers in the labels correspond to 
the numbers on atoms on the left panel. As for the eigenvalue, we only show for the range 
[—0.01,0.01] with color scale shown on the right and the contours for 0.01 and —0.01 are 
shown by white dashed lines. The triangle dots shows the locations of the Lagrange points. 
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Figure 5: The largest eigenvalue of the stress tensor and corresponding eigenvector of 
AI4H10, plotted in the same manner as Fig. [4j 
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Figure 6: The largest eigenvalue of the stress tensor and corresponding eigenvector of 
AI4H12, plotted in the same manner as Fig. [4j 
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AI4H10 





Figure 7: The largest eigenvalue of the stress tensor and corresponding eigenvector of 
AI4H10 (but on the different plane from Fig. [51), plotted in the same manner as Fig. [il 
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